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Abstract — Owing  to  the  favorable  ambiguity  function  properties 
and  the  increased  deployment,  mobile  communications  systems 
are  useful  for  passive  bistatic  radar  applications.  Further,  simulta¬ 
neously  using  multiple  illuminators  in  a  multistatic  configuration 
will  improve  the  radar  performance,  providing  spatial  diversity 
and  increased  resolution.  We  compute  modified  Cramer-Rao 
lower  bounds  (MCRLB)  for  the  target  parameter  (delay,  Doppler) 
estimation  error  using  universal  mobile  telecommunications 
system  (UMTS)  signals  as  illuminators  of  opportunity  for  passive 
multistatic  radar  systems.  We  consider  both  coherent  and  non-co- 
herent  processing  modes.  These  expressions  for  MCRLB  are  an 
important  performance  metric  in  that  they  enable  the  selection  of 
the  optimal  illuminators  for  estimation. 

Index  Terms — Coherent  processing,  Cramer-Rao  bound,  dis¬ 
tributed,  multistatic,  passive  radar,  UMTS  signals. 


I.  Introduction 

PASSIVE  radar  systems  use  several  signals  of  opportunity 
as  illuminators  for  target  estimation  unlike  their  active 
counterparts  that  require  expensive  transmission  equipment. 
Some  examples  of  these  signals  include  television  [1],  [2], 
audio  broadcast  signals,  FM  radio  [3],  [4]  and  mobile  commu¬ 
nications  systems  [5].  By  employing  just  the  receivers,  these 
systems  are  inherently  robust  and  not  detectable  [6]— [8].  With 
their  wide  spread  availability,  mobile  communications  signals 
are  important  illuminators  of  opportunity.  Further,  employing 
passive  radar  systems  in  a  multistatic  configuration  provides 
spatial  diversity  and  improves  the  resolution  [9]  similar  to 
the  concept  of  active  MIMO  radar  systems  [10]— [13].  These 
passive  multistatic  radar  systems  view  the  targets  from  different 
aspect  angles  illuminated  by  different  transmitters. 

UMTS  is  a  third  generation  wireless  communications  stan¬ 
dard.  These  signals  have  been  shown  to  provide  good  resolu- 
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tion  as  an  illuminator  for  passive  bistatic  radar  systems.  [14] 
covers  a  detailed  ambiguity  function  analysis  using  these  illu¬ 
minators  in  bistatic  configuration  to  study  the  global  resolution 
performance  in  terms  of  mainlobe  width  and  sidelobes.  Further, 
it  also  computes  the  bistatic  MCRLB  for  these  systems.  CRLB  is 
important  for  analyzing  the  local  estimation  accuracy  as  it  pro¬ 
vides  a  lower  bound  on  the  error  variance  of  unbiased  estima¬ 
tors.  Further,  under  certain  conditions,  the  maximum  likelihood 
estimators  (MLE)  asymptotically  achieve  the  CRLB  [15]— [17]. 
While  non-coherently  combining  the  data  from  different  con¬ 
stituent  bistatic  pairs  provides  spatial  diversity  gain,  whenever 
phase  synchronization  is  possible,  these  systems  can  be  used  for 
coherent  processing  to  obtain  very  high  resolution  properties. 

In  this  paper,  we  will  compute  the  MCRLB  [1 8]  on  the  target 
estimation/localization  using  UMTS  based  passive  multistatic 
radar  under  both  coherent  and  non-coherent  processing  sce¬ 
narios.  Since  the  transmitted  symbols  are  not  deterministic,  it  is 
difficult  to  compute  the  classical  CRLB.  Therefore,  we  shall  use 
the  MCRLB  as  an  alternative  approach  here  by  averaging  the 
Fisher  information  matrix  over  the  probability  mass  function  of 
the  transmitted  symbol.  In  [19],  [20],  MCRLB  has  been  studied 
and  applied  to  radar  problems.  In  [14],  the  authors  computed  the 
MCRLB  for  passive  bistatic  radar  using  UMTS  signals.  For  the 
non-coherent  processing  mode,  we  will  compute  the  modified 
Fisher  information  matrix  (MFIM)  by  incorporating  the  UMTS 
passive  bistatic  radar  results  of  [14]  into  the  CRLB  analysis  for 
non-coherent  MIMO  radar  described  in  [21].  For  the  coherent 
processing  mode,  in  this  paper  we  will  derive  closed-form 
expressions  for  the  MFIM  by  computing  the  expected  values  of 
all  the  second  order  derivatives  of  the  log-likelihood  function. 

Our  MCRLB  analysis  will  provide  a  quantitative  measure 
of  passive  multistatic  radar  performance.  The  expressions  will 
not  just  be  a  function  of  the  transmitted  waveforms  but  also  a 
function  of  the  multistatic  radar-target  geometry.  The  geometry 
has  been  shown  to  play  a  very  important  role  in  determining 
the  bistatic  ambiguity  function  [22],  thereby  also  impacting  the 
bistatic  CRLB.  While  passive  radar  does  not  offer  the  freedom 
of  designing  the  transmitted  waveforms,  it  does  allow  flexibility 
in  selecting  the  transmitters  from  among  several  choices.  The 
geometry-dependent  MCRLB  analysis  will  open  up  a  new  di¬ 
mension  for  passive  multistatic  radar  systems  by  aiding  the  se¬ 
lection  of  optimal  illuminators  of  opportunity  to  achieve  a  de¬ 
sired  target  estimation  accuracy. 

The  rest  of  the  paper  is  organized  as  follows.  In  Section  II, 
we  present  the  non-coherent  and  coherent  measurement  model 
for  passive  multistatic  radar.  In  Section  III,  we  will  compute 
the  expressions  for  the  MFIM  for  the  non-coherent  scenario  fol¬ 
lowed  by  the  coherent  scenario  in  Section  IV.  We  will  present 
some  numerical  examples  to  demonstrate  our  analytical  results 
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in  Section  V.  We  shall  demonstrate  that  the  MCRLB  is  signif¬ 
icantly  lower  for  the  coherent  processing  case  than  the  non-co- 
herent  processing  case.  Further,  we  will  demonstrate  the  depen¬ 
dence  of  the  MCRLB  values  on  the  multistatic  geometry.  Fi¬ 
nally,  we  provide  concluding  remarks  with  potential  future  work 
in  Section  VI. 


II.  Signal  Model 

Consider  a  system  comprising  of  MT  transmitters  and  Mr 
receivers.  Let  the  ith  transmitter  be  located  at  ti  =  \txi,tyi]. 
Similarly,  the  jth  receiver  is  located  at  rj  =  [rxj,ryj].  The 
target  location  and  velocity  are  represented  by 

~P=\Px,Py],  (1) 

'v>=[vx,vy\.  (2) 

Note  that  we  chose  the  vectors  from  the  2-dimensional  Cartesian 
space  for  simplicity.  The  results  can  easily  be  extended  to  the 
3 -dimensional  space  without  loss  of  generality.  Therefore,  we 
define  the  target  state  vector 

M  =  \Px,Py,vx,vy]  .  (3) 

Note  that  in  the  rest  of  the  paper,  we  restrict  our  analysis  to 
the  single  target  scenario.  In  the  presence  of  multiple  targets, 
the  number  of  unknown  variables  in  the  target  parameter  vector 
increases  by  a  factor  equal  to  the  number  of  targets,  thereby 
making  the  analysis  much  more  complex.  In  future  work,  we 
will  extend  the  analysis  to  the  multiple  target  scenario  by  ap¬ 
pending  the  parameters  corresponding  to  multiple  targets  into 
the  parameter  vector.  The  baseband  signal  corresponding  to  the 
ith  transmitter 


1  ^ 

Ui(t)  -  -y=  ^2  cin9i{t  -  nT),  (4) 

n=0 

where  c;n  are  the  transmitted  quadrature  phase  shift  keying 
(QPSK)  symbols,  N  is  the  number  of  symbols  and  T  is  the 
inverse  of  the  chip  rate.  The  pulse  gi{t)  is  defined  as  delayed 
root-raised-cosine  (RRC)  pulse  gi(t)  =  hi  (t  —  y),  where  the 
delay  D  and  h{(t)  are  defined  in  [14]. 

Let  r\lrj  and  denote  the  different  bistatic  delays  and 
Doppler  shifts  associated  with  the  target  located  at  fi : 


where  c  represents  the  speed  of  wave  propagation  in  the  medium 
and  fc  represents  the  carrier  frequency.  These  transformation 
relations  between  the  delay-Doppler  space  and  the  Cartesian  co¬ 
ordinates  are  necessary  while  computing  the  MCRLB. 


angle  of  view  and  the  widely  spaced  antennas  decorrelate  these 
attenuations  into  uncorrelated  zero-mean  complex  Gaussian 
random  variables.  If  the  target  RCS  values  corresponding  to 
certain  bistatic  pairs  are  weak,  it  is  highly  likely  that  they  will 
be  compensated  for  by  other  pairs  with  strong  returns.  We  will 
recollect  the  expressions  for  the  multistatic  radar  measurement 
vector  using  non-coherent  processing.  Note  that  the  signals 
from  different  transmitters  are  assumed  to  be  separable  at  the 
receivers  in  some  domain  (for  example  different  frequency 
spectra).  Therefore,  we  have  MtMr  components  of  the  re¬ 
ceived  signals. 


Vijit)  =  PijUi  (t  -  Tij)  e 


3  2tt 


+  (7) 


Note  that  y^t)  are  statistically  independent  for  different  trans¬ 
mitter-receiver  pairs  because  of  the  wide  separation  between 
the  antennas  that  leads  to  independent  looks  of  the  target.  Let 
riij  (t)  denote  the  additive  zero-mean  white  Gaussian  noise  with 
variance  a2 .  Further,  the  target  attenuations  faj  are  zero-mean 
Gaussian  distributed  with  variance  a2. 

Note  that  the  transmitted  symbol  is  not  deterministic.  There¬ 
fore,  we  will  average  the  Fisher  information  matrix  over  the 
probability  mass  function  of  the  transmitted  symbol  c.  First,  for 
a  given  c,  the  log-likelihood  ratio  across  all  transceiver  pairs 
[21] — [24] 


ATt  -V/r 

log/(yV)|c)  =  (s/iiWiq ,  (8) 

i—  1  j—  1 


where 


1  (VSMIC)  = 


/2=3,o-w“**(t-ro-)e 


xe 


d  t 


Therefore, 


loS^  |c)  = 


r2  O2  + 


Xui  (*  -  r!j) e 


/oo 

vm 

-oo 


d  t 


+  C, 


where  C  is  independent  of  the  target  state  vector  fi.  Note  that 
the  carrier  dependent  phase  term  is  not  present  because  it  is  ab¬ 
sorbed  into  the  circularly  symmetric  Gaussian  distributed  atten¬ 
uation  terms  and  does  not  impact  the  non-coherent  processing. 

Given  this  statistical  model,  the  expression  for  the  Fisher  in¬ 
formation  matrix  for  non-coherent  MIMO  radar  with  widely 
separated  antennas  was  derived  in  [21].  Using  the  relationships 
between  the  delay  terms  and  the  Cartesian  coordinates, 


III.  Non-Coherent  MCRLB 

In  this  section,  we  will  compute  the  MCRLB  for  the  non-co¬ 
herent  processing  scenario  by  deriving  the  MFIM  expression. 
Under  this  scenario,  the  target  is  made  up  of  several  individual 
isotropic  scatterers  [11].  The  target  attenuations  vary  with  the 


^^~ij  _  1  |  Px  ^xi  ^  Px  Gc j  j 

dpx  C  \j  _  t*\\  || "p  -  rj\\)  ’ 

f  Py^  ~  ^ yi  _|_  Py  ~  ry:i  ) 

dpy  ~c  VllV-TlI  IIG-Ullj' 
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Similarly,  the  relations  between  the  Doppler  shift  terms  and  the 
Cartesian  coordinates  for  target  position  and  velocity  give  [2 1  ] 


and 


d/g, 

dpx 


~vx 

A 
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- - - b 

||V  -  t,  ||  ||p  -  Tj  | 

Px 


df£. 

dpy 


A||“p  -  U  || 

{vx  (txi  -  px)  +  Vy  ( tyi  -  Py )) 

j  ^  xj  Px 

'  ~  ,  3 

All  p  -  rill 

(vx  (rxj  -  px )  +  vy  (ryj  -  py)) , 
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- — - b  _ ^  ^ 

||V  -  U  II  II V  -  rj | 

tyi  ~  Py 


A||~p  -  Hi  || 

(Wx  ( txi  -  Px)  +  Vy  ( tyi  -  Py)) 

■  ryj  ~  Py 

'  _  _ ^  3 

M\  p  -  rj\\ 

(vx  {rxj  -  px)  +  Vy  (ryj  -  Py))  , 


where  A  denotes  the  carrier  wavelength,  and 


dvx 

df^ 
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^  xj  Px 


A||  p"  - 


dvy 


txi  Px 

Any  -  ~tl 

tyi  ~  Py 

Any -^ii  Aiiy-^ 
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Py 


The  analysis  for  non-coherent  MIMO  radar  in  [2 1  ]  expresses 
the  MIMO  FIM  as  a  combination  of  the  constituent  bistatic 
FIMs.  Therefore,  we  can  express  the  MFIM  for  our  problem  as 


J(M)  = 


87 r2cr4 


ATx  A/r 


(m): 


(9) 


where  the  final  expressions  for  the  elements  of  the  symmetric 
Fisher  information  matrix  Jij  (fi)  corresponding  to  the  ijth 
transceiver  pair  are  given  in  Appendix  using  the  bistatic  results 
from  [14]  and  the  equations  for  the  transformation  of  variables 
described  above.  Note  that  since  the  transmitted  waveforms 
U{(t)  are  not  deterministic,  expected  values  of  the  entries  of 
FIM  were  used  to  arrive  at  the  expressions  in  the  Appendix. 
This  is  in  contrast  to  the  classical  CRLB  that  needs  us  to 
average  across  the  transmitted  waveform  space  in  the  joint 
probability  density  function  of  the  measurement  vector  y^(t) 
instead  of  the  conditional  density  function.  Computing  the 
classical  CRLB  is  not  feasible  for  our  problem  [14].  However, 
MCRLB  has  been  studied  and  applied  to  radar  problems  earlier 
[19],  [20]  as  it  provides  a  good  benchmark  when  computing 
the  classical  CRLB  is  not  feasible.  Finally  we  will  invert  this 
modified  FIM  to  obtain  the  modified  CRLB  as  the  diagonal 
entries  of  J  1  (/1).  We  clearly  observe  from  the  expressions 
of  the  entries  of  the  modified  FIM  (in  the  Appendix)  that 
the  MCRLB  depends  on  the  multistatic  geometry  along  with 
the  transmitted  waveform  parameters  ai ,  thereby  making  the 
choice  of  transmitters  important. 


IV.  COHERENT  MCRLB 

Having  computed  the  MCRLB  for  the  non-coherent  pro¬ 
cessing  scenario,  we  move  on  to  the  coherent  processing  mode. 
In  this  mode,  the  target  attenuations  are  made  up  of  a  single 
point  scatterer  that  has  an  isotropic  complex  reflectivity  denoted 
by  (3  =  PRe  +  jPi m[10].  Further,  this  attenuation  coefficient  is 
a  deterministic  unknown  quantity  unlike  the  random  Gaussian 
modeling  in  the  non-coherent  case.  Therefore,  we  modify 
the  target  state  vector  to  include  the  attenuation  coefficients 
(i  —  [px,py,vx,vy,  /3 Re,  Am]-  An  alternative  representation 
of  the  target  state  vector  in  the  delay-Doppler  domain  is 
y—  rfl5  .  .  .  ,T^-tMr,  /^ii5  .  .  .  ,  Re,  Am  .  We  will 

use  the  representation  given  by  v  as  an  interim  step  before 
computing  the  MCRLB. 

Since  we  have  assumed  that  the  different  transmitted  signals 
are  separable  at  the  receivers  (different  carrier  frequencies/re¬ 
sources),  the  received  signal  corresponding  to  the  ith  transmitter 
and  the  jth  receiver  for  the  coherent  processing  scenario  can  be 
given  as 

Vij(t)  =  Put  (t  -  rft)  e  V  ,J  '3)  +ni:i(t). 

Note  the  presence  of  the  phase  term  e_j27r^cir"  in  the  above  ex- 
pression.  The  CRLB  for  coherent  MIMO  radar  with  stationary 
targets  was  computed  in  [25]-[27].  Recently,  moving  targets 
were  considered  in  [28]. 

The  log-likelihood  function  across  all  receivers  for  a  given 
transmitted  waveform 

AEp  A/r 

iogi  (yv)ic)  =  y^yyog  1  (^(f)ic) ,  (io) 

i=  1  .7  =  1 

where 


-  72  JZc  Vij  «  ('-  U  )  ' 


d  t 


Therefore, 
log  1  (y^it)  |c) 

Mt  Mr  «oc 
n  i=  1  j=  1 


This  simplifies  to 

log  l  (y^(t)\c)  = 

Mt  Mr  rOQ 


Vi  jit)  -  Pui  {t  -  rtj) 

J27r  (/d  .  .  it~Tij  )—fciTij  ) 


dt  +  C. 


1 


lvlit  poo  / 

H  (  iviM* Vijit) 

i= 1  j= 1  v 

*  (t  —  t?.  )  —  /ciTM  ^ 

,  ,^,2  *  /.  /,  M 

+  \P\  ui  [t  ~  Tij)  ui  [t  ~  Tij)  e  v  3  ' 

(/£  .^  (*-^ )  ) 

,  .  -jilt  (fit  (t-r^)-fc  irA\ 

-FtfjitKit-Tfte  V-VA  -7  Mjdi+c, 

where  C  is  a  constant  independent  of  the  target  state  vector  (i. 
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The  third  term  reduces  to 

MT  poo 

\Pf  /  u*  (; t  -  rfi)  Ui  [t  -  t£)  d t  =  |/3|2.  (11) 

i= 1 

Note  that  we  are  considering  unit  energy  complex  envelopes 
for  the  QPSK  transmitted  waveforms  shaped  by  RRC  filters  as 
given  in  [14].  Further,  the  first  term  is  purely  a  function  of  the 
measurements  and  does  not  contain  the  delay,  Doppler,  and  at¬ 
tenuation  parameters.  Hence, 


Mt  Mr  poo  ,  , 

logl(yIJ(t)\c)  =  ( ReV 

n  i= 1  j= 1  J~oo  \  l 


Ui  (*  ~  Tij)  e 

mtmr\p\ 2 


d  t 


+  c, 


where  C'  is  another  constant. 

It  is  easier  to  compute  the  MFIM  in  the  v  domain  and  then 
apply  a  change  of  variables  transformation  to  the  /i  space. 
Therefore,  we  have 


J(M)  = 


(12) 


First,  we  need  to  compute  the  term  that  corresponds  to  the 
change  of  variables.  In  this  term,  we  have  already  stated  the  ex¬ 
pressions  for  the  derivatives  of  the  delay  and  Doppler  terms  with 
respect  to  the  Cartesian  positions  and  velocities  in  the  non-co- 
herent  processing  section.  The  terms  that  are  new  for  this  co¬ 
herent  scenario  are  the  derivatives  with  respect  to  the  attenua¬ 
tions  Prg  and  Pim.  Since  the  delays  and  Dopplers  do  not  depend 
on  these  terms,  these  derivatives  are 


=  0. 


K 


dp  Re  dpi 

Additionally,  we  also  have 


=  0, 


dJk 

dpRe 


=  0, 


dpi  rr 


=  0. 


dPRe  _  q  dPim  _  o  dPRe  _  o  dPim 
dpx  ’  dpx  7  dpy  7  dpy 

dPRe  _  Q  dpi  m  _  0  d Prq  _  0  dpi m 

dvx  7  dvx  7  dvy  7  dvy 

dpRe  1  dpRe  dPim  0  9  Pi  m 

dPRe  ~  '  9Pim  ~  7  dPRe  ~  ’  dPim 


Now,  the  only  matrix  we  need  to  compute  for  obtaining  the 
MCRLB  is  J{v).  We  computed  closed- form  expressions  for 
each  entry  of  this  matrix  by  evaluating  all  the  second-order 
derivatives(see  Appendix  for  the  derivations).  Using  these  ex¬ 
pressions,  we  will  compute  the  coherent  MCRLB  as  the  diag¬ 
onal  elements  of  MCRLBc  (//)  =  J1  (/jl). 

V.  Numerical  Examples 
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Fig.  1 .  Simulated  multistatic  scenario  with  Transmitter  and  receiver  locations. 


system.  We  consider  5  transmitters  and  5  receivers  (See  Fig.  1) 
located  at 


t\  —  [1.5,  2]  km  ;  rf  =  [3,0]  km; 

T?  =  [1.5,  2.25]  km  ;  rj  =  [5,0]  km; 

=  [1.5,  2.5]  km  ;  rj  =  [0,3.5]  km; 

~tl  =  [2, 1.5]  km  ;  rt  =  [1, 1]  km; 

=  [2, 1.75]  km  ;  F5  =  [4,  4]  km. 

We  will  compute  the  square  root  of  MCRLB  (RMCRLB)  around 
the  position  [3,  4]  km  and  velocity  [30,  50]  m/s.  We  chose  the 
same  system  parameters  as  in  [14]  for  the  simulations;  obser¬ 
vation  time  NT  =  0.1s,  T  =  0.26  /is,  a  =  0.22,  and  the 
center  frequency  fc  =  2100  MHz.  Define  the  signal-to-noise 
ratio  (SNR)  as 

SNR  =  10  log  (13) 

In  Fig.  2,  we  plotted  the  non-coherent  RMCRLB  in  the  x-po- 
sition  and  y-position  dimensions  as  a  function  of  the  SNR.  Sim¬ 
ilarly,  Fig.  3  shows  the  velocity  RMCRLB  for  varying  SNR.  We 
observe  that  the  RMCRLB  is  lower  in  the  y-dimension  for  both 
the  position  and  velocity.  At  an  SNR  of  0  dB,  the  RMCRLB  for 
the  x  and  y  positions  are  9.6206  m  and  4.8496  m,  respectively. 
The  RMCRLB  for  the  x  and  y  velocities  at  the  same  SNR  are 
0.1727  m/s  and  0.1121  m/s,  respectively.  Now,  to  show  the 
importance  of  the  geometry,  we  change  the  position  for  which 
we  are  computing  the  RMCRLB  to  [3, 1.5]  km.  While  this  does 
not  affect  the  terms  E  {ei},  E  {7^},  and  E  {r]ij},  it  impacts 
the  derivatives  of  the  delay-Doppler  terms  with  respect  to  the 
Cartesian  coordinates.  We  clearly  observe  from  Fig.  4  that  the 
RMCRLB  at  0  dB  SNR  for  px  is  5.5077  m  and  that  for  py  is 
5.1567  m.  These  numbers  are  different  from  the  earlier  case  and 

the  same  holds  true  even  for  the  velocity  RMCRLB.  _ 

For  the  coherent  processing  scenario,  we  choose  p  =  . 

Therefore,  the  signal-to-noise  ratio 


In  this  section,  we  will  use  numerical  examples  to  compute 
the  MCRLB  for  a  UMTS-based  passive  multistatic  radar 


SNR  =  10  log 


(14) 
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Fig.  2.  Non-coherent  RMCRLB  in  the  x-position  and  y-position  dimensions  as 
a  function  of  SNR. 


Fig.  3.  Non-coherent  RMCRLB  in  the  x- velocity  and  y-velocity  dimensions  as 
a  function  of  SNR. 


Fig.  4.  Non-coherent  RMCRLB  in  the  x-position  and  y-position  dimensions  as 
a  function  of  SNR  when  [px,py]  =  [3,1.5]  km. 


We  plotted  the  RMCRLB  curves  using  the  same  parameters  as 
the  non-coherent  case  around  the  position  [3,  4]  km  and  velocity 


Fig.  5.  Coherent  RMCRLB  in  the  x-position  and  y-position  dimensions  as  a 
function  of  SNR. 


Fig.  6.  Coherent  RMCRLB  in  the  x-velocity  and  y-velocity  dimensions  as  a 
function  of  SNR. 


[30,  50]  m/s.  In  Fig.  5,  we  plot  the  RMCRLB  for  both  the  po¬ 
sition  parameters.  We  clearly  notice  that  the  RMCRLB  is  much 
lower  when  compared  with  the  non-coherent  case. 

Next,  we  plotted  the  RMCRLB  for  the  velocity  parameters 
in  Fig.  6.  Even  here,  the  CRLB  is  significantly  lower  when 
compared  with  the  non-coherent  scenario.  At  an  SNR  of 
—  6  dB,  the  coherent  velocity  RMCRLB  in  x  and  y  dimensions 
are  0.3076  m/s  and  0.1237  m/s,  respectively  as  opposed  to 
0.5460  m/s  and  0.3545  m/s  for  the  non-coherent  case.  For  the 
coherent  case,  there  are  additional  unknown  parameters  other 
than  the  positions  and  velocities  in  the  form  of  /3.  In  Fig.  7,  we 
plot  the  RMCRLB  for  /JRe  and  (3im  as  a  function  of  the  SNR. 
Just  as  we  did  for  the  non-coherent  case,  we  will  change  the 
position  to  different  value  and  demonstrate  the  dependence  of 
the  RMCRLB  values  on  the  geometry.  Let  the  true  position 
parameters  be  [6,  2.5]  km.  Fig.  8  demonstrates  the  change  in 
values  of  RMCRLB  when  compared  to  the  previous  geometry 
in  Fig.  5.  When  the  true  position  is  [6,  2.5]  km,  we  observe 
that  the  x-position  dimension  has  a  lower  RMCRLB  than  the 
y-position  dimension.  This  is  contrary  to  the  scenario  when  the 
true  target  position  was  [3,  4]  km.  In  [14],  the  authors  show  that 
the  bistatic  RMCRLB  shoots  up  rapidly  as  the  target  position 
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have  considered  both  the  non-coherent  and  coherent  processing 
scenarios.  We  observe  that  the  MCRLB  is  significantly  lower  for 
the  coherent  processing  case  than  the  non-coherent  processing. 
Further,  we  demonstrated  the  dependence  of  the  MCRLB  values 
on  the  multistatic  geometry.  The  MCRLB  results  validate  the 
feasibility  of  using  the  UMTS  signals  as  an  important  illumi¬ 
nator  for  passive  radar. 

In  future  work,  we  will  extend  this  research  by  consid¬ 
ering  other  illuminators  of  opportunity  in  addition  to  UMTS 
waveforms.  Additionally,  we  will  extend  the  analysis  to  the 
multiple  target  scenario  by  appending  the  parameters  corre¬ 
sponding  to  multiple  targets  into  the  parameter  vector  and 
evaluating  all  the  corresponding  additional  terms  in  the  MFIM. 
Further,  this  would  facilitate  a  way  to  solve  the  important 
problem  of  optimal  illuminator  selection  for  passive  multistatic 
radar  systems. 


Appendix 


Fig.  8.  Coherent  RMCRLB  in  the  x-position  and  y-position  dimensions  as  a 
function  of  SNR  when  [px,py]  =  [6,2.5]  km. 


approaches  the  bistatic  baseline.  In  the  multistatic  case,  several 
constituent  bistatic  pairs  affect  the  overall  RMCRLB  and  hence 
it  is  more  complex  to  study  the  dependence  of  these  expressions 
on  the  geometry. 

Note  that  even  though  the  coherent  processing  mode  offers 
lower  RMCRLB  than  the  non-coherent  mode,  these  two  sce¬ 
narios  represent  two  completely  different  target  models.  There¬ 
fore,  coherent  processing  is  not  applicable  when  the  target  con¬ 
sists  of  multiple  individual  isotropic  scatterers  that  cause  uncor¬ 
related  phase  shifts  across  different  transceiver  pairs.  Further, 
it  is  very  important  to  have  perfect  phase  synchronization  be¬ 
tween  all  the  transmitters  and  receivers  for  coherent  processing 
and  this  is  not  always  possible  as  a  result  of  several  physical 
limitations  like  inaccurate  knowledge  of  the  antenna  locations 
and  local  oscillator  characteristics  [29]— [3 1  ]. 

VI.  Concluding  Remarks 

We  have  computed  MCRLB  for  passive  multistatic  radar  sys¬ 
tems  using  the  UMTS  waveforms  as  the  illuminators  of  oppor¬ 
tunity  by  deriving  closed-form  expressions  of  the  MFIM.  We 


The  entries  of  the  constituent  bistatic  MFIM  for  the  non-co¬ 
herent  processing  scenario  can  be  expressed  as 


Jilin) 

^(m) 
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JiiM 
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For  the  coherent  processing  scenario,  first,  we  will  compute 
aiosy)lc) .  Recollect  from  Section  IV  that  the  log-likelihood 

ij 

log!  (^(0 1  c) 


MtMrW  |  c< 

/t-2 


We  observe  that  the  first  and  third  terms  in  the  above  expression 
will  be  zero  as  they  do  not  depend  on  the  delay  terms.  Hence, 


Next,  we  compute  the  derivative  with  respect  to  the  Doppler 
shift  term.  Since  the  Doppler  term  only  shows  up  in  the  carrier 
and  not  in  the  waveform  term,  we  have 


Finally,  the  first  derivative  w.r.t  the  attenuation  terms 
dlog  l  <yy)  |c)  _  _  2 MTMR/jRe 

dpRe  <*n 


and 

dlogl  (y^(t)\c)  _  _  2MTMRf3im 

dpim 


We  clearly  observe  that 

d2  logCjCMM 

^Re2 
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Further, 

<92  log l  (y^(t)\c)  47 r 
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To  compute  the  MFIM,  we  need  to  take  the  expected  value 
of  these  second  order  derivative  terms  using  the  joint  density 
function  of  the  measurements  and  the  transmitted  codewords. 
Therefore, 

(  d2  logl  (yv(t)\c)\ 
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Since  the  transmitted  symbols  at  different  intervals  are  indepen¬ 
dent,  we  obtain 
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From  the  first  term,  we  notice  that  when  compared  to  the 
second  term,  it  has  the  variable  t  inside  the  integral.  Therefore, 
this  term  becomes  an  odd  function  and  integrates  to  zero.  The 
second  term 
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We  have  finished  computing  first  order  derivatives.  However, 
the  FIM  needs  the  evaluation  of  the  second-order  derivatives. 
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Hence, 
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After  simplification 

52log/(2/'i(/)|c)l  __ 


E 


dfLdHiu 


-E<% 
a: 


poo 

J  I  Im  (i  -  7%)  (t  -  t£)  u*  (t  -t£)}  d t 

n  J  —oo 


Using  similar  process  as  the  results  above,  we  obtain 
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Therefore,  taking  the  expected  value  of  this  expression,  we  get 
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It  can  be  easily  shown  that  this  equation  simplifies  to 
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The  first  term  in  the  above  expression  is  zero  because  the  trans¬ 
mitted  RRC  pulses  are  even  functions  and  the  first  derivatives 
of  even  functions  are  always  odd  functions,  thereby  forcing  the 
value  of  the  two  sided  integral  to  zero.  Therefore, 
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Taking  the  expected  value  of  this  term 
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Following  the  same  approach  as  earlier,  the  first  term  of  the 
above  expression  becomes  zero  because  the  term  inside  the  inte¬ 
gral  is  an  odd  function.  Further,  the  waveforms  from  each  trans¬ 
mitter  are  unit  energy  waveforms.  Therefore,  we  have 
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So  far,  we  have  computed  all  the  second  order  derivatives  in¬ 
volving  the  attenuations  variables.  The  second  order  derivatives 
with  respect  to  the  Doppler  terms 
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The  term  inside  the  integral  in  the  above  equation  has  been 
derived  in  [14]  while  computing  the  MCRLB  for  the  passive 
UMTS-based  monostatic  radar  systems. 


E 


{/: 


t2  I Ui  {t)\2  dt  ^  = 


T2  D2 

16 al  +  T~ 

DT(N  —  1)  T2(N  -  1)(2N  -  1) 


Substituting  this  result,  we  finally  express  the  second-order 
derivatives  with  respect  to  the  Doppler  terms  as 
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Next,  we  have  the  second-order  cross  derivative  terms  in¬ 
volving  both  the  delay  and  the  Doppler  frequency  parameters 


Upon  substituting  the  expressions  for  the  expected  values  of  the 
measurements,  we  obtain, 
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We  can  simplify  the  above  expression  to 
d2  log  l  \c)  \  _ 


+ 


8*  (/.I  +  /£,,)  r  (  f 


E 


i:n 


ij\  ;  ft  v 

ij 


Kdfo, 

poc  N  —  1 


/°°  /  (  / 

E  Re  U'  (t  +  nT  + 

-oo  ^—n  V  l  V 


47 r 

N(Jn  r)  -oo  n=Q 


x  e 


"j27r  (^,,(*-^5)-^^) 


dt 


dh*  ( t )  1  \ 

x  '■•w  a  L,d 

/‘oc  A^  —  1 


8tt2  (/ci  +  /£,.) 


EH 


P*Vij(.t)ui  {t  -  Tij) 


A  pOO  iy  —  ± 

+  JTZ 2  /  E  (Re  {■?’  l^l2  ^  (*)  h*  0)})  di 

n  J  - 00  n=0 

rNfuj 

n  — n  A  ^ 


x  e 


dt. 


Ncl 

x  hi  (t)  h*  (i)  jj  di 
This  further  simplifies  to 
<92  log  l  (yA‘(0|c) 


n=0 


l/E  ( t  +  nT  +  — 


The  expected  value 

d2  log  l(y^(t)\c) 


E  ■ 


4 E  (  ^  Re  { \fS\2  m  (t  -  t'E)  TA  J,  d t 

® n  I  •/  oc 


E 


dTtjdfk 


87 r 


47 r 


iV_1  /  r  rc 

■E  Re  ££2  / 


thi  (t) 


d_K  (0 

dt 


+ 


4ir 


} 

P— -L  (  /*  OO 

E  Re  y  i£2  /  hi  (o  hi  wdi  [ 

n=0  ^  >'-00  J 

ra=0  '  ^ 

J  hi  (t)  h*  (i)d/j^ 


+ 


(/ci  +  /£ ) 


*3 


r( 


Re<i|/3|2Mi  (t-rfi) 


,  du*  (*  -  r£) 


*3 


^  ^  '  dt 


^  n=0 

87T2  ( fc\  +  / n..  )  ^Z.1  /  f  /  £) 

n<t2  ^l"~ll/31  (nT  +  V 

n=0 


87T 


+  -- 


(/d  +  //£ ) 


/  (Re  { l^l2  Wi  (*  -  Tij) 

J  —  00  v  ^ 


X<  (i-£')})di}- 


The  second  term  in  the  above  expression  contains  an  odd  func¬ 
tion  in  the  integral.  Therefore, 


Note  that  we  used  the  fact  that  the  first  derivative  of  the  RRC 
waveform  is  an  odd  function.  Further,  for  the  term,  the  expres¬ 
sion  inside  the  integral  is  purely  imaginary  because  the  RRC 
waveform  is  real.  Therefore, 
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where  Hi{uj)  denotes  the  Fourier  transform  of  hi(t).  Here,  we 
used  the  result  in  equation  (53)  of  [14]  for  obtaining  the  last 
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